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I. INTRODUCTION 

The nucleon-nucleon (NN) interaction problem has a 
well deserved reputation of being a difficult one (for a 
review see e.g. |lj). One obvious reason may be found 
in the theoretical and experimental inaccessibility of the 
shortest distance interaction region relevant to Nuclear 
Physics. This lack of fundamental knowledge might be 
remedied in the future by ah initio lattice calculations, 
for which some incipient results exist already 0, Q . How- 
ever, the current situation does not prevent from address- 
ing many facets of NN interactions, particularly long dis- 
tance features, provided short distance insensitivity is 
guaranteed. Despite the great efforts during the years 
based on successful phenomenology [1, IE H, 0] , only dur- 
ing the last decade has the parentage to the underlying 
QCD been made more explicit after the proposal [8|, |9( 
to use Effective Field Theory (EFT) methods and ex- 
tensive use of Chiral S y mmetr y (for comprehensive re- 
views see e.g. Ref. [M El El)- EFT starts out with 
the most general Lagrangean in terms of the relevant 
degrees of freedom compatible with chiral symmetry, di- 
mensional analysis and perturbative renormalization, or- 
ganized by a prescribed power counting based on assum- 
ing a large scale suppression on the parameters 4^/^ ~ 
A/jv ~ lGeV. Thus, by construction, EFT complies to 
the expectation of short distance insensitivity, provided 
that enough counter-terms encoding the unknown short 
distance physics are added in perturbation theory to a 
finite order 0, EE Q EE EE El EE EE HJ ■ However, 
the physics of bound states, such as the deuteron, is gen- 
uinely non-perturbative and an infinite resummation of 
diagrams proves necessary (see however [2l| for a renewed 



perturbative set up). There arises the problem of non- 
perturbative renormalization and/or modifications of the 
original power counting for which no universally accepted 
scheme has been found yet, mainly because the require- 
ments may vary. Thus, although the EFT method does 
comply to model independence it is far from trivial to 
achieve regulator independence and a converging pattern 
dictated by a reasonable power counting within a non- 
perturbative set up. 

The NN renormalization problem may be cast ef- 
ficiently in the traditional language of potentials and 
the corresponding Schrddinger equation. Besides phras- 
ing the NN interaction into the familiar and more in- 
tuitive non-relativistic quantum mechanical framework, 
this procedure has also the advantage that the pertur- 
bative determination of the potential complies to the 
desired short distance insensitivity for the potential it- 
self @, llQEQil. The unconventional feature 
is that chiral expansions necessarily involve singular po- 



tentials at short distances, i.e. r \V(r)\ 



for 



0. 



Indeed, in the limit r«l/ (or equivalently large mo- 
menta) pion mass effects are irrelevant in loop integrals 
and hence at some fixed order of the expansion one has 
for the coupled channel potential 1 



V(r 



M N c 2n 



+m+r 



(47r/ 7r ) 2n M™A r r 2n + m + r 



(1) 



where n, m and r are nonnegative (> 0) integers, A the 
NA splitting, and Cfe, with k = In + m + r, is a dimen- 
sionless matrix of van der Waals coefficients in coupled 
channel space. The dimensional argument is reproduced 
by loop calculations in the Wein berg dimensional power 
counting @, E3, E3, EE EE 123, l23f l23| and is scheme 
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1 The only exception is the singlet channel-One-Pion-Exchange 
(OPE) case which behaves as ~ m^/ f%r, see below. 
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independent. It is thus conceivable that much of our un- 
derstanding on the physics deduced from chiral potentials 
might be related to a proper interpretation of these highly 
singular potentials, which degree of singularity increases 
with the order of the expansion. This of course raises 
the question about in what sense higher order poten- 
tials are smaller. One obvious way to ensure this desired 
smallness is by keeping a finite and sufficiently moderate 
cut-off d,[2l,|2l,[ll[2i,[13,[li|, so the singularity of the 
potential is not probed, effectively recovering the power 
counting expectations for the size of higher order con- 
tributions. The disadvantage is that results are strongly 
cut-off dependent for scales about r c ~ 0.5 — 1.0 fm, sim- 
ilar to the ones probed in NN scattering (H, [3(J HJ H2] • 
However, this is not the only possibility. Explicit com- 
putations [29|, [3(| [3l| show that, when the cut-off is re- 
moved, the higher order contributions will continue to 
generate only small changes in the amplitudes under spe- 
cific circumstances which can be determined a priori. As 
emphasized in previous works [2ij , [30L [3~fl | renormaliza- 
tion is the most natural tool provided 1) we expect the 
potential to be realistic at long distances and 2) we want 
short distance details to be inessential in the descrip- 
tion. As discussed above, this is precisely the situation 
we face most often in Nuclear Physics. A surprising and 
intriguing feature is that knowledge on the attractive or 
repulsive character of the singularity, i.e. the sign of the 
eigenvalues of in Eq. |T]), turns out to be crucial to suc- 
cessfully achieve this program and ultimately depends on 
the particular scheme or power counting used to compute 
the potential. 

The singularities of chiral potentials may be discon- 
certing 2 , but they can be handled in a way that do not 
differ much from the standard treatment of well-behaved 
regular potentials which one usually encounters in nu- 
clear physics [35j]. Renormalization is the mathematical 
implementation of the appealing physical requirement of 
short distance insensitivity and hence a convenient tool 
to search for typical long distance model- and regulator 
independent results. In a non-perturbative setup such 
as the NN problem, renormalization imposes rather tight 
constraints on the interplay between the unknown short 
distance physics and t he p erturbatively computable long 
distance interactions [29l l30l . l3ll ] . This viewpoint pro- 
vides useful insights and it is within such a framework 
that we envisage a systematic and model independent 
description of the NN force based on chiral interactions. 
In this regard it is amazing to note how the sophisticated 
machinery of perturbative renormalization in Quantum 
Field Theory used to compute the chiral potentials has 
not been so extensively developed when the inevitable 
non-perturbative physics must be incorporated; quite of- 
ten the renormalization process is implemented by trial 



2 Early treatments can be found in |33j (see Ref. [34| for an early 
review) . 



and error methods by adding counter-terms suggested 
by the a priori power counting in momentum space. 
Howe ver, detailed analyses in coordinate and momentum 
space [li, US HIM Ha, HI] show that the allowed struc- 
ture of counter-terms can be anticipated on purely ana- 
lytical grounds and cannot be chosen independently on 
the long distance potential (see also Refs. [3§, H(| for nu- 
merical work). Specifically, for a channel-subspace with 
good total angular momentum the number of countert- 
erms is n(n + l)/2 with n the number of negative eigen- 
values of the van der Waals matrix Cfc in Eq. (JTJ) . In fact, 
the much simpler coordinate space renormalization has 
been shown to be fully equivalent to the popular momen- 
tum space treatment for purely contact theories WlJ and 
theories containing additional long range physics [421 |43| . 
On the light of these latter studies the smallness of in- 
creasingly singular potentials such as Eq. ((T|) is trig- 
gered quite naturally by choosing the regular solution 
of the Schrodinger equation; the wave function behaves 
as u p (r) ~ (^Kf^r) ? (modulo prefactors depending 
on the attractive or repulsive character of the potential). 
Thus, as shown in Ref. [431 ]. when the short distance cut- 
off r c approaches a fixed scale, ~ 1/(4^/^) = 0.2fm, an 
increasing 0( r ™ +m / 2+r//2+1 ) insensitivity of phase shifts 
and deuteron properties is guaranteed as the power of 
the singularity of the potential increases. Indeed, calcu- 
lations with chiral Two Pion Exchange (TPE) potentials 
reproducing low energy NN data display this insensitivity 
for reasonable scales of r c ~ 0.5fm [29|, [3(| which corre- 
spond to the shortest wavelength probed by NN scatter- 
ing in the elastic region. Within such a scenario, the 
discussion on whether one should remove the cut-off or 
not [1^, [3(J 0, HU would become less relevant as 
the order of the chiral expansion is increased in the com- 
putations. This requires of course that the same renor- 
malization conditions are implemented, in order for the 
computations to represent an specific physical situation, 
and that the cut-off lies in the stability region r c ~ 0.5fm, 
which in turn means that renormalization has effectively 
been achieved. 

However, despite all these findings, the question on 
how a sensible hierarchy for NN interactions should be 
organized has been left open. Although we know whether 
and, in positive case, how this can be made compatible 
with the desired short distance insensitivity [lij y(| [H| , 
not all chiral potentials based on any given power count- 
ing are necessarily eligible. The Weinberg counting based 
in a heavy baryon approach at leading order (LO) [ll| 
for 1 Sq and 3 Si — 3 Di states turns out to be renor- 
malizable (29j due to an attractive-repulsive character of 
the coupled channel eigen potentials at short distances. 
There is at present no logical need why this ought to 
be so, for the simple reason that power counting does 
not anticipate the sign of the interaction at short dis- 
tances. In fact when one goes to next-to-leading order 
(NLO) the short distance 1/r 5 singular repulsive char- 
acter of the potential makes the deuteron unbound (30j 
because the interaction becomes singular and repulsive. 
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Finally, next-to-next-to-leading order (NNLO) potentials 
diverge as — 1/r 6 and are, again, compatible with Wein- 
berg counting in the deuteron [30], in this case because 
the interaction is singular and attractive. Further incon- 
sistencies between Weinberg's power coun ting and renor- 
malization have been reported in Refs. [32, E2]. The 
previous examples show clearly that the requirement of 
renormalizability can be in open conflict with the idea 
of a convergent pattern based on a preconceived chiral 
expansion and the mere dimensional power counting of 
the original proposal [8J. In particular, the assumption 
that the NN potential cannot be completely determined 
at arbitrary short distances is incompatible with a fully 
repulsive short distance singularity because finitcness im- 
plies that no counter-term is allowed. 3 

In this regard one should note an interesting analogy 
between the NN interaction in the chiral quark model 
(for a review see e.g. (4(| H3, |48j]) and the van der Waals 
molecular interactions in the Born-Oppcnhcimcr approx- 
imation [3(| • For non-relativistic constituent quarks the 
NN interaction is provided by the convoluted OPE quark- 
quark potential in the direct channel since quark ex- 
change and finite nucleon size effects can be neglected 
for distances larger than the nucleon size. Second or- 
der perturbation theory in OPE among quarks generates 
TPE between nucleons yielding 

V NN = (NN\V OP e\NN) 

+ y \(NN\V OPE \HH>)\ 2 | (2) 

hJ^nn Enn ~ Ehh ' 

for the NN potential, where |iVJV) represents a two- 
nucleon state and \HH') an arbitrary intermediate state. 
In Regge theory one has the relation M\ — M N = 
m 2 — m 2 [49J, m p being the rho mass, which in the chi- 
ral limit and for Mjy ~ N c m p /2 implies a NA split- 
ting A - m 2 p /2M N ~ m p /N c = 256MeV for N c = 3 
and suggests a scale numerically comparable to twice 
the pion mass, which actually vanishes for N c — > 00. 
Thus, it makes sense to consider an EFT where the AN 
splitting is regarded as a small parameter [H(| HH • Ac- 
tually, when HH' = NA and HH' = AA, Eq. © 
resembles the result found using the Feynman graph 
technique [Til ]. Moreover, the second order perturba- 
tive character suggests that the potential becomes at- 
tractive, since Enn — Ena — —A and Enn — Baa = 
— 2A. At short distances the matrix elements scale as 
(NN\V pe\HH') ~ g 2 A /(f 2 r 3 ) and hence the potential 
becomes singular ~ —g\/(Af^r 6 ) and attractive, neces- 
sarily being renormalizable with an arbitrary number of 
counter-terms through energy dependent boundary con- 
ditions [36| . Clearly, the renormalization of a potential 
where the NA splitting is treated as a small scale de- 



3 In the coupled channel case that means all eigen-channels being 
repulsive at short distances. 



serves further investigation since it appears as an obvious 
candidate where all necessary requirements for a conver- 
gent and short distance insensitive scheme might be met. 
In the present paper we check that the naive expecta- 
tions based on the simple Eq. ([2]) are indeed verified for 
the long distance chiral potentials including A degrees of 
freedom [|, [H Q H3 . 

The delta-isobar has played a crucial role in the devel- 
opment of nuclear and particle physics (see e.g. [H, HU 
and references therein). Besides fixing the number of 
colours N c — 3 in QCD to comply with Pauli princi- 
ple both at the quark as well as at the hadron level, 
this state can be clearly seen in irN scattering as a res- 
onance, and it ubiquitously appears whenever the nu- 
cleon excitation energy is about the Delta-Nucleon split- 
ting, A = Ma - M n = 293MeV. Already in the ear- 
liest implementations of Weinberg's ideas the influence 
of Delta degrees of freedom was considered [1, [H, [56| 
in old-fashioned perturbation theory where energy de- 
pendent potentials are generated. Energy independent 
potentials have been obtained using the Feynman dia- 
gram technique a decade ago [HI, [lj] and only recently 
the N2LO contributions have also been worked out in 
Ref. [521 ] . where peripheral np phase shifts are computed 
in perturbation theory for these A contributions. Fi- 
nite cut-off calculations involving A-degrees of freedom 
to NLO have been analyzed in momentum space (25|. It 
has been shown that the inclusion of the A-excitation im- 
proves the con verg ence of the chiral expansion of the NN 
interaction [57ll58|. The renormalization of the 1 Sq phase 
with the NLO-A chiral potential has also been discussed 
recently in momentum space (42[ . Short distance insensi- 
tivity has also been discussed in other contexts different 
from NN collisions such as nd scattering at threshold [59[ 
where a rationale for the multiple scattering expansion 
was indeed supported by the chiral singular potentials. 
Unfortunately, boost corrections being proportional to 
the average relative pn kinetic energy in the deuteron, 
{p 2 )d/M, turned out to diverge due to the infinitely many 
short oscillations of the wave function. Recently, the ben- 
efits of including the A in such a process are discussed 
in [60| where it is shown how a proper reorganization of 
the boost corrections not only makes them finite but also 
numerically small after renormalization in harmony with 
phenomenological expectations. The drawback is a pro- 
liferation of low energy constants when including the A, 
which can render the A-full theory impractical at higher 
orders. The implications of a A-based power countin g fo r 
the three nucleon problem are analyzed in Refs. [57l l58j|. 

In the present paper we focus our interest in the crucial 
role played by the Delta-isobar in NN scattering in the 
elastic region from the point of view of renormalization 
of chiral nuclear forces, and how it might solve a long- 
standing problem. Besides an acceptable phenomenology 
in the s-wave phases we also show that precisely because 
of the built-in short distance insensitivity and unlike pre- 
vious calculations where the A-degrees of freedom were 
absent, the deuteron is always bound up to N2LO-A, the 
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highest order computed at present [521 ]. 

The paper is organized as follows. In Section [II] we re- 
view the formalism as applied to the energy independent 
potentials in a A-full theory [HQ HI- In Section ITTT1 
we discuss the simpler 1 Sq channel including either one 
or two counter-terms by means of an energy dependent 
boundary condition. We analyze the deuteron bound 
state and its properties in Section llVl Scattering states in 
the 3 Si — 3 Di channel as well as the corresponding phase 
shifts are constructed by orthogonality to the deuteron 
and compared to the results obtained with the Nijmegen 
II potential @, which has a \ 2 P er datum near one and 
therefore can be considered as an alternative partial wave 
analysis, compatible with the original Nijmegen PWA 5]. 
The results found in this paper are further compared to 
the A-less theory in Section [V] with previous results [3C| . 
In Sec. |VI] we address the issue of the relevant scales en- 
tering the calculation as well as the role of missing effects 
such as relativistic, 37r contributions or vector meson ex- 
change calculations. We also re-analyze a regularization 
of the potential based on the s pec tral representation mo- 
tivated by dispersion relations [26|, |27|, |28|, |52j and discuss 
its meaning on the light of the present approach. Finally 
in Section I VIII we recapitulate our results and draw our 
main conclusions. 



II. FORMALISM 

In the present paper we use coordinate space renormal- 
ization by means of boundary conditions [23. l30l l3ll] . The 
equivalence to momentum space renormalization using 
counter-terms for s ing ular potentials is discussed thor- 
oughly in Refs. [13, E3| . We use the e nerg y inde p endent 
long distance chiral potentials of Ref. [lj, [lj, [52| , which 
are on-shell equivalent to the energy-dependent poten- 
tials of Ref. We follow the convention of taking 1/M 
corrections as 0(p 2 ), in agreement with the original ar- 
gument from Ref. [6l| as well as with the calculations of 
Refs. H, HI, and the computation of the N2LO-A po- 
tential of Ref. [52]. We remind that this convention is 
not universal, and keeping the 1/M corrections as 0{p) 
is also customary 4 . This convention has the consequence 
that at NLO-A and N2LO-A relativistic corrections are 
suppressed [52j |. The final expressions for the potential 
are taken from Ref. [Hj]. 

We solve the coupled channel (coupled in angular mo- 
mentum) Schrodinger equation for the relative motion, 
which in compact notation reads, 



M 



+ V NN (x) 



(3) 



4 Although the effect of keeping or ignoring the mass corrections 
is indeed small (see end of this section). 



where the spin and isospin indices have not been explic- 
itly written; M = 2M p M n /(M p + M n ) is twice the re- 
duced proton-neutron mass. In coordinate space the po- 
tential can be written as 



V NN (x) = Vc(r)+TW c (r) 

+ a{V s {r)+TW s {r)) 
+ S 12 (V T {r)+TW T (r)) 



(4) 



where spin-orbit and quadratic spin-orbit terms have 
been ignored as they are not present in the NLO-A and 
N2LO-A potentials from Ref. [52T |. The operators r, a 
and S12 are given by 



t = n-T 2 = 2t(t + i) - 3, 

a = d\ ■ c?2 = 2s(s + 1) — 3 , 
5i2 = 3 cti • fa 2 • f — 3\ ■ <?2 , 



(5) 



where fi(2) and c?i(2) represent the proton(neutron) 
isospin and spin operators; t\ -Tz and <j\ -a 2 are evaluated 
for total isospin t = 0, 1 and total spin s — 0, 1. Note 
that in Ref [52j the potential is local and there are no 
relativistic mass corrections at N2LO-A. For states with 
good total angular momentum j and total spin s = 1, 
the tensor operator reads 



/ 2(j-l) 

f)» 11 



°12 



2j+l > 





6^(3 + 1) 
2j+l 



- 



2(3+2) 
2j+l 



(6) 



where the matrix indices represent the orbital angular 
momentum / = j — 1, j, j + 1. For s = states the tensor 
operator vanishes S{ 2 = 0. We remind that Fermi-Dirac 
statistics implies (— ) l+s+t = —1. 

A remarkable feature which happens at LO, NLO-A 
and N2LO-A @, [H, [3 [13] is that the spin-orbit cou- 
pling vanishes, as well as any relativistic corrections. As 
a consequence, for a given isospin the potential can be di- 
agonalized. The corresponding eigen-potentials depend 
just on the isospin of the channel and not on the total 
angular momentum j. All the j dependence goes into 
the matrix which diagonalizes the potential. This can be 
seen in the following formulas: 

V y = (V c + tW c ) - 3 (V s + W s t) + S{ 2 (V T + tW t ) 
= Mj diag (A -4B,A + 2B, A + 2B) Mj x , (7) 

where V lj is the triplet-channel potential written in ma- 
trix form, with 

A = (Vc + tWc)-3(Vs + W s t), 
B = {Vt + tWt), 

^cos 9j — sin 8j \ 
Mj = 1 , (8) 
ysin^- cos Qj J 
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W riNLO-A 

IVIU 6,W,S 




u/iNLO-A 
IvlLj 6,W,T 


Fit-1 


-7.233 


-1.201 


0.601 


0.301 


-0.601 


-0.301 


Fit-2 


-3.867 


-0.833 


0.417 


0.177 


-0.417 


-0.177 




a / r/-tN2LO-A 
1V1L, 6,V,C 


a / r/-tN2LO-A 
1V1L, 6,W,C 




aj- ,-<N2LO-A 
IVIL/ 6,W,S 




i /fr <N2LO-A 
JWL/ 6,M/,T 


Fit-1 


-1.241 


-0.489 


0.244 


0.351 


-0.244 


-0.351 


Fit-2 


-4.381 


-1.121 


0.561 


0.439 


-0.561 


-0.439 



TABLE I: van der Waals MC& coefficients (in fm 4 ) for the different spin-isospin components of the NLO-A and N2LO-A 
potentials. We use the nN motivated Fits 1 and 2 of Ref. [&|. Fit 1 involves the SU(4) quark- model relation Ka = 3gA/(2y/2) = 
1.34 for g A = 1.26. 



where 

cosg 7 = a J . (9) 

V 2j + 1 

The transformation Mj diagonalizes the full potential but 
not the Schrodinger equation, since it contains in addi- 
tion to the potential the centrifugal barrier, which is a 
diagonal operator in the jls basis, but does not remain 
diagonal in the rotated basis. Specific knowledge of the 
short distance behaviour is needed to carry on with the 
renormalization program. Generally, on purely dimen- 
sional grounds, we have for r — > 

f-iV,i s~tW,i 

Wi(r)--£-, (10) 

where 

fj* \ nil 

L/ fc=2n+m+r+l j2n ]y[ m £±r ' \ I 

with A the iVA splitting, and n, m and r nonnegative 
integers. At short distances, the angular momentum de- 
pendence may be neglected when the index k > 2. The 
relevant issue to carry out the renormalization procedure 
and to generate finite results is to know whether the in- 
teraction is attractive or repulsive. It turns out that both 
NLO-A and N2LO-A potentials have a leading singular- 
ity behaviour of 1/r 6 . In Appendix [Al we list the analyt- 
ical expressions for the van der Waals coefficients CI , for 
k = 6 (i.e. the leading singularity of the potential). 

In our numerical calculations we take / w = 92.4MeV, 
m ff = 138.03MeV, 2^ np = M N = 2M p M n /(M p + M n ) = 
938.918MeV, g A = 1-29 in the OPE piece to account 
for the Goldberger-Treimann discrepancy and 17,4 = 1.26 
in the TPE piece of the potential 5 . The correspond- 
ing pion nucleon coupling constant takes then the value 



5 Strictly speaking gA = 1-26 both in the OPE and TPE pieces 
of the potential, but it happens that at NLO and higher orders 
the OPE piece receives a contribution from the dig LEC, related 
to the Goldberger-Treimann discrepancy, see, for example, the 
expressions of Ref. [26| for details. This is equivalent to consider 
the original expression for the OPE potential, but taking gA = 
1.29 instead of g A = 1-26. 



9ttNN — 13.083 for the OPE piece of the potential. We 
use the AN splitting A = 293 MeV. For h A and the low 
energy constants Ci , C2 , C3 , C4, 63 and 6g we take the val- 
ues from Fit 1 and 2 of Ref. [52| (tabic I within that ref- 
erence) , where they are deduced from a fit to irN thresh- 
old parameters in S- and P- waves to the data of Ref. [13] • 
These values are compatible with all ttN threshold pa- 
rameters except 6^ _j_ which is about twice its recom- 
mended value [62| at this level of approximation. Fit 
1 uses the SU(4) quark-model relation Ka = 3gA/2V% 
(= 1.34 for g A = 1.26), while Fit 2 uses h A = 1.05. 

The numerical values of the van der Waals coefficients 
MCq are summarized in Table [J for the different compo- 
nents of the NLO-A and N2LO-A potentials. The van 
der Waals coefficients are additive: therefore one can ob- 
tain the coefficient corresponding to some given partial 
wave by adding the individual contributions, i.e. 

MC 6 = MC 6 , v ,c(r)+TMC 6tW ,c 
+ o-(MC 6 ,v, s + tMCg,w,s) 
+ S 12 (MC 6 ,v,t + tMC 6 ,w,t) ■ (12) 

This is done for the singlet 1 5o and the triplet 3 Si — 3 Di 
channels in Table HT1 

We also present N2LO-^ results for comparison pur- 
poses. For them we use the same parameters as in the 
N2LO-A computation, except for the c\, C3 and C4 LECs, 
for which we use the following two different determina- 
tions: the so-called set IV of Ref. [30| , which was ob- 
tained in Ref. [24[ b y fi tting to NN scattering data, and 
the values from Ref. [52| for N2LO-^ which we refer to as 
set 7rN and allow a better comparison with the N2LO-A 
computations presented here, as they are also fitted to 
reproduce ttA^ S- and P-wave threshold parameters. Due 
to the different fitting procedures, any direct compari- 
son between the results of Ref. [30], i.e. set IV, and the 
N2LO-A results should be done with care. It should be 
mentioned too that the results of Ref. [3(| contain 1 /M 
corrections to the potential. These corrections are nev- 
ertheless small and, if excluded, would only induce small 
differences in the results of Ref. |f0l] 6 . 



As an illustration, ignoring the 1/M corrections will yield to 
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6/Sn 










-Rt 


Fit 1 


-11.138 


-3.932 


0.855 


-4.537 


-5.141 


-3.327 


Fit 2 


-6.481 


-1.482 


0.322 


-1.710 


-1.938 


-1.254 




«,^N2LO-A 
M °6,K5n 




MC 6 N |°- A 


MC R N f n °- A 


-i4 


-Rt 


Fit-1 


-3.519 


-0.585 


2.290 


-2.204 


-3.823 


1.035 


Fit-2 


-8.500 


-1.773 


2.138 


-3.284 


-4.796 


-0.261 



TABLE II: van der Waals MC 6 coefficients (in fm 4 ) in the and 3 Si - 3 Di (deuteron) channels of the NLO-A and N2LO-A 
potentials. We also present the corresponding negative eigenvalues —R\ and — R*l in the deuteron channel case. We use the 
?riV motivated Fits 1 and 2 of Ref. JH]. Fit 1 involves the SU(4) quark-model relation h A = 3£u/(2V2) = 1.34 for g A = 1.26. 



III. THE SINGLET CHANNEL 



where 



A. Equations and boundary conditions 

The 1 Sq wave function in the pn center-of-mass (cm.) 
system can be written as 



47T r 



(13) 



with the total spin s — and m s = 0. The function 
u(r) is the reduced S- wave function, and satisfies the 
following reduced Schrodinger equation 

-u'l(r)+U 1So (r)u k (r) = k 2 u k (r) , (14) 

where k is the center-of-mass momentum, and Uis the 
reduced potential defined as 



CTi So (r) = M (V c (r) + W c (r] 
-3V s (r) - 3W 8 {r) 



(15) 



For a finite energy scattering state we solve for the chiral 
potential with the asymptotic normalization 



Uk(r) 



sin(fcr -I- So(k)) 
sin<5o(fc) 



(16) 



with 8o(k) the phase shift. For a potential falling off 
exponentially ~ e - m ^ r a t large distances, we have the 
effective range expansion, valid for momenta |fc| < 



kcot5 Q (k) = -— +-r k 2 +v 2 k 4 
ao 2 



(17) 



with ao the scattering length and the effective range. 
At short distances the NN chiral potential behaves as 



MC 6 i So 



R* 

r 6 



(18) 



MC, 



e^So 



MC; 



NLO-A 



MC; 



N2LO-A 



(19) 



which is a van der Waals type interaction 7 . The nu- 
merical values for MC K N i?" A and MC^ff~ A are listed 
in Tab. HU The value of the coefficient is negative, with 
the typical length scale R = (—MCq) 1 / 4 . The solution 
at short distances is of oscillatory type: 



Uk{r) 



, , r \3/2 . 
A[-) sm 



1 (R 



2 7' + ^ 



(20) 



where A is a normalization constant and tp an undeter- 
mined phase, which may in principle depend on energy. 

We will present below two calculations. In the first 
one the renormalization is carried out with one counter- 
term, which in turn means one renormalization condition. 
In such a case the short distance phase becomes energy 
independent and orthogonality conditions are satisfied. 
In the second calculation we proceed with two counter- 
terms, i.e. two renormalization conditions, for which the 
short distance phase acquires a very specific energy de- 
pendence. 



B. Renormalization with one counter-term 

As mentioned, the phase shift is determined from 
Eq. (|16p . but to fix the undetermined phase ip we im- 
pose orthogonality for r > r c between the zero energy 
state and the state with momentum p. As shown in [301 ] , 
orthogonality turns out to be equivalent to the follow- 
ing condition between 0- and k-momentum reduced wave 
functions at r = r r 



u' k (r c )u (r c ) - u' (r c )u k (r c ) = . 



(21) 



Taking the limit r c — ► implies that the short distance 
phase (p is energy independent [30l ]. Thus, for the zero 
energy state we solve 



<(r) + Ui So (r)u (r) = 0. 



(22) 



the following modifications for the results of Ref. 13011 (origi- 
nal results in parentheses): A$ = 0.887(0.884) fm 1 / 2 , r m = 
1.972(1.967) fm, Q d = 0.278(0.276) fm 2 , P D = 8(8)%, (r" 1 } = 
0.442(0.447) fm -1 , <r" 2 ) = 0.276(0.284) fm" 2 . 



7 It should be noted that the NLO-A. potential is less singular than 
the NLO-A one at short distances, diverging as 1/r 5 only. 
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k cm [MeV] (c) k cm [MeV] (d) 



FIG. 1: Renormalized phase shifts for the LO, NLO-A and N2LO-A potentials as a function of the cm. np momentum fc c . m . 
in the 1 So singlet channel, compared to the Nijmegen II potential results The computations are done both with Iia = 1.34 
(upper panels) and with Ka = 1-05 (lower panels). For the LO result, the calculation is always done with one counterterm, 
while for the NLO and N2LO results, computations are done with one counter-term (left panels) and two counter-terms (right 
panels), corresponding this last case to the standard Weinberg counting for the 1 So channel. 



with the asymptotic normalization at large distances 

Mr) -»!-—, (23) 



where ao is the scattering length. In this equation ao 
is an input, so one needs to integrate Eq. (|22p from in- 
finity to the origin (in contrast with the usual procedure 
of integrating from the origin to infinity). The effective 
range, defined as 



To 



2 / dr 
lo 



1 



r 



u {rf 



(24) 



can be computed. Due to the superposition principle, 
we can write the zero momentum wave function as the 
following linear combination 



u (r) = uo, c (r) - — «o,«(r) 
ao 



(25) 



where uo, c (r) — * 1 and uo_ s (r) — ► r correspond to cases 
where the scattering length is either infinity or zero re- 
spectively. Using this decomposition, we get 



r = A 
where A, B and C, defined as 
.4 



B 



C 



a a\ ' 



2 / dr(l - ^ ) , 
Jo 

POO 

B = -4 / dr(r - u 0tC u 0lS ) , 
Jo 



C 



dr(r 2 - ul ) . 



(26) 

(27) 
(28) 
(29) 



depend on the potential parameters only. The interesting 
thing is that all dependence on the scattering length ao 
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TABLE III: Predicted threshold parameters in the singlet 1 So 
channel with one counter-term for Fits 1 and 2 of Ref. [52] • We 
compare our renormalized results given by the cut-off inde- 
pendent universal formula (|26p for ro and its extension for V2 
to finite cut-off NN calculations using their scattering length 
as an input. The experimental values for the scattering length 
and effective range are taken from Ref. [7J. 





Calculation 


ao(tm) 


ro(rm) 


V2 (im ) 


LO 


Ref. [30] 


Input 


1.44 


-2.11 


NLO-A 


Ref. [30] 


Input 


2.29 


-1.02 


NLO-A {h A = 1.34) 


This work 


Input 


2.91 


-0.32 


NLO-A (h A = 1.05) 


This work 


Input 


2.76 


-0.53 


N2LO-A (Set IV) 


Ref. [30] 


Input 


2.87 


-0.38 


N2LO-A (ttN) 


This work 


Input 


2.92 


-0.31 


N2LO-A (Fit 1) 


This work 


Input 


3.05 


-0.12 


N2LO-A (Fit 2) 


This work 


Input 


3.04 


-0.13 


Nijm II 


Ref. [5, 6] 


-23.73 


2.67 


-0.48 


Reid 93 


Ref. [5, 6] 


-23.74 


2.75 


-0.49 


Exp. 




-23.74(2) 


2.77(5) 





is displayed explicitly by Eq. 
r = 2.661- 
r = 2.517- 
r = 2.789- 
r = 2.780- 



5.707 


5.988 


a 


al 


5.430 
! 


5.811 


a 


al 


5.992 
1 


6.187 


OiQ 


a 2 


5.969 
1 


6.171 


a 





Numerically we get 
(NLO-A, h A = 1.34) , 
(NLO-A, h A = 1.05) , 
(N2LO— A, Fitl) , 
(N2LO— A, Fit2) . 

(30) 



The corresponding numerical values when the experimen- 
tal ao = — 23.74fm is taken, as well as the V2 parame- 
ter can be looked up in Table IIIII As can be seen, the 
value of the effective range has a clear convergence pat- 
tern when going from LO to N2LO-A. The contribu- 
tion coming from N2LO-A is very small compared to the 
NLO-A contribution, in agreement with the findings of 
Ref. [52j. This is in contrast with the A-less theory, for 
which N2LO generates a great correction over the NLO 
results. Unfortunately the N2LO-A results converge to 
wrong values, about 3fm for both Fit 1 and 2. The same 
happens for the N2LO-^ results, although with a weaker 
convergence pattern. The reason for this discrepancy is 
that the NLO and N2LO potentials, both in A-less and 
A-full theories, are too attractive at intermediate dis- 
tances, thus yielding a bigger value for tq than the one 
obtained with phenomenological potentials like Nijmegen 
II or Reid93 HIH. 

Using the orthogonality condition, Eq. ([21]) . the phase- 
shift can be determined from the scattering length and 
the potential as independent parameters when the limit 
r c — > is taken. The renormalized phase shift is pre- 



sented in Fig. [T] (left). The phase shifts are compared 
with the ones obtained from the Nijmegen II potential 0, 
which are compatible with the Nijmegen PWA As we 
see the trend in the effective range tq and the vi param- 
eter is reflected in the behavior of the phase shift. In 
Ref. [42| a similar calculation was carried out in momen- 
tum space with inclusion of A degrees of freedom at NLO 
and one counter-term. The present NLO-A coordinate 
space results agree with that calculation when the SU(4) 
relation, tiA = 1-34, is used. As discussed in the previous 
paragraph for the effective range parameters, at N2LO- 
A the phase shifts have already converged, although to 
the wrong value, due to the excessive strength of the in- 
termediate range of the potential. In the next section we 
will see how the situation changes when an extra counter- 
term is included in the computations. This counter-term 
will be fitted to reproduce the effective range. 



C. Renormalization with two counter-terms 

In the standard Weinberg counting both the long dis- 
tance potential as well as the short distance potential are 
dimensionally expanded in momentum space. For the 
short range piece, one writes Vs(p',p) — Co + C 2 {p 2 + 
p 12 ) + . . . , where Co, C2, etc, are referred to as counter- 
terms. In Ref. [42j the interrelation between the counter- 
terms and short distance boundary conditions in coor- 
dinate space has been discussed at length. At LO in 
the Weinberg counting one has only one counter term. 
This corresponds to the situation described in the previ- 
ous section. At NLO and N2LO in the Weinberg count- 
ing (both in the A-less and A-full cases), one adds two 
counter-terms Co and C2 for a fixed momentum space 
cut-off A, which may be fixed by reproducing the scat- 
tering length ao and the effective range ro . The interest- 
ing finding was that such a procedure does not generate 
a converging 1 S'o phase shift in the limit A — > 00 
Thus, the momentum space polynomial parameterization 
of Weinberg counting is not compatible with renormal- 
ization. However, this does not necessarily mean that 
one cannot impose two renormalization conditions to fix 
the scattering length and the effective range. Actually, 
on a wider perspective one may pose on the one hand the 
problem of obtaining a finite phase shift embodying the 
chiral potential and on the other the problem of fixing 
both ao and ro as independent parameters. Fortunately, 
there exists a unique procedure in coordinate space [36| 
meeting the two previous conditions and yielding conver- 
gent phase shifts, which we apply below to discuss the 
^o channel when two counter-terms are considered for 
NLO-A and N2LO-A chiral potentials. We will refer to 
this scheme as Weinberg counting with A. 

The equivalent coordinate space procedure [H, |42j 
consists of expanding the wave function in powers of the 
energy 



Uk(r) = Uo(r) + k 2 u 2 {r) 



(31) 
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where ito(r) and U2(r) satisfy the following equations, 
-uff(r) + tfi So (r)uo(r) = 0, (32) 

/ \ i — >oo T 

u (r) 1 , 

-u>>(r) + U 1So (r)u 2 (r) = u Q (r) , (33) 

a (r 3 - 3a r 2 + 3a r r) 



u 2 (r) 



6a 



The asymptotic conditions correspond to fixing a and 
ro as independent parameters (two counter-terms). The 
matching condition at the boundary r = r c becomes en- 
ergy dependent [36[ 



u' k (r c ) _ u' (r c ) + k 2 u' 2 (r c ) + 



Uk(r c ) u Q (r c ) + k 2 u 2 (r c ) + 



(34) 



whence the corresponding phase shift may be deduced 
by integrating in Eq. (j3"2"|) and Eq. (|33|) and integrat- 
ing out Eq. ([TJ} with Eq. ([To]) . It is worth mentioning 
that the energy dependent matching condition, Eq. (fM)) . 
is quite unique since this is the only representation for 
the boundary condition guaranteein g th e existence of the 
limit r c — > for singular potentials [30|. As pointed out 
in Refs. [36, 42], polynomial expansions in k 2 such as sug- 
gested e.g. in the Nijmegen PWA and implemented 
later on for chiral TPE potentials do not work for 
r c — ► 0, and in fact generate undesired oscillations for 
r c -c 1.4fm in the phase-shifts. The validity of these fea- 
tures can be deduced analytically in coordinate space as 
a consequence of the RG-invariance of a Moebius bilinear 
transformation [36j | . Equivalent parameterizations in mo- 
mentum space may likely exist, but are so far unknown. 
As already mentioned, the widely used polynomial rep- 
resentation of short distance interactions in momentum 
space V s (p',p) = C + C 2 (p 2 +p' 2 ) + . ■ ■ of standard NLO 
and NNLO Weinberg counting implies that for A — > oo 
either C 2 is irrelevant when only ao is kept fixed or the 
phase shift does not converge when both ao and ro are 
simultaneously fixed [H, HH • 

In Fig. [1] (right panel) we show the results for the two 
counter-term renormalized phase shift at NLO-A and 
N2LO-A. As we see, the second counter-term is respon- 
sible for a certain improvement all over the elastic region, 
in particular in the region k < where TPE effects are 
expected to dominate in the singlet 1 So -channel 8 . 

At higher momenta, however, the discrepancy with the 
Nijmegen II potential results [y] and therefore with the 
Nijmegen PWA fB| persists. The second counter-term 
is therefore unable to provide the necessary repulsion to 



In principle for k < m-jr OPE (instead of TPE) effects should 
dominate, but in the case of the 1 Sq singlet channel this is not 
the case due to the weakness of OPE in this channel. This only 
happens in the singlet channel; in the case of the triplet channel 
we recover the naive expectations, and OPE clearly dominates 
for k < . 



compensate for the excessive intermediate range attrac- 
tion present in the NLO-A and N2LO- A potentials. This 
result agrees with the findings of Ref. [42j and extends 
them from NLO-A to N2LO-A. Possible solutions to this 
disturbing situation include finite cut-off computations, 
which help to increase the effect of the second counter- 
term, and the inclusion of spectral function regulariza- 
tion, proposed in Ref. [26| precisely to treat the problem 
of the intermediate strength of the N2LO-A potential in 
peripheral partial waves. The case of spectral function 
regularization will be discussed in Section IVII 



IV. THE TRIPLET CHANNEL 

A. Equations and boundary conditions 

The 3 Si — 3 Di wave function in the pn cm. system 
can be written as 

1 



w(r) 



u(r)a p ■ a n 



V8 



(3cr p ■ xa n ■ x — a p ■ a n ) 



X 



pi I 



(35) 



with total spin s = 1 and m s — 0, ±1; a p and <j n are 
the Pauli matrices for the proton and the neutron re- 
spectively. The functions u(r) and w(r) are the reduced 
S- and D-wave components of the relative wave function 
respectively. They satisfy the coupled set of equations in 
the 3 Si — 3 Di channel 

-u"{r) + U3 Sl {r)u[r)+U El (r)w(r) = k 2 u{r) , 



-w"(r) + U El (r)u(r) 



£W»0 + 4 



w(r) = k 2 w(r) , 
(36) 

with U3 Sl (r), Ue 1 ( , t) and Uau x {f) the corresponding ma- 
trix elements of the coupled channel potential, which are 

U 3Sl = V c - -iW c + V s - 3W S , 
U El = 2V2(V T -3W T ), 
Us Dl = v c -3Wc + Vs-3Ws-2Vt + 6W t . 



(37) 



At short distances one has 



Us 



MC, 



Si 



6. 3 Si 



r 

MC, 



6 

6, Ei 



r 

MC, 



6 ' 



U El 

Us Dl 



where the van der Waals coefficients are given by 
Cqs Si 



(38) 



/-.NLO | ,-<N2LO 
°6, 3 Si + °6, 3 Si i 



n _ /-(NLO , /-fN2LO 

^6,£i — <^6, Ei + °6,Ei 7 



c 



6, 3 Di 



/-.NLO | /-.N2LO 
6, 3 E! + °6, 3 Di ■ 



(39) 
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FIG. 2: Deuteron wave functions, u (left panel) and w (right panel), as a function of the distance (in fm) for the LO, NLO-A 
and N2LO-A potentials, compared to the Nijmegen II wave functions [6J. The asymptotic normalization u(r) — + e~ lr has been 
adopted and the asymptotic D/S ratio is taken r\ — 0.0256(4). The oscillations of the wave functions are related to the presence 
of deeply bound states, but do not have any appreciable effect on deuteron observables as they cannot be resolved in effective 
field theory. 



Their numerical values are listed in Table |TT] for Fits 1 
and 2 of Ref . [52j . One can diagonalize the corresponding 
matrix of van der Waals coefficients 

(mc 6 3 Si mc 6>Ei \ 

\MC 6 , El MC 6j3D J 



cos t) 


sin 


— sin ( 


1 cos 


cos 9 


— sin 


sin 9 


cos 6 




(40) 



where, according to Eq. ([£ 
cos 9 



the angle is 
1 

: 7s' 



(41) 



i.e. 9 = 54.7° and common to both the NLO-A and 
N2LO— A potentials. That means that the eigenvalues 
of the NLO-A and N2LO-A van der Waals matrices are 
additive, and therefore can be summed up directly from 
Table HU 



(42) 



L ± — Jl -±,NLO Jl ±,N2LO ■ 

As we see from Table HH the NLO-A matrix is negative 
definite, but N2LO-A matrix is not so for Fit 1. How- 
ever, what counts is the NLO-A + N2LO-A contribution 
which, as can be checked, is negative definite. In the di- 
agonal basis one has at short distances 




cos 9 sin 
— sin 9 cos 

where the short distance eigen functions are 

2 r 2 



(43) 



v + (r) 



r 



C+ sin 



v-(r) = ( — ) C-sin 



2T 2 " 



<P- 



(44) 



and ip± are short distance phases which must be fixed 
independently on the potential and C± suitable normal- 
ization constants. Orthogonality of solutions of different 
energy requires these phases to be energy independent. 
Following the procedure of Ref. [30[ we fix them from 
deuteron physical properties, namely the binding energy 
and asymptotic D/S ratio (see below). Once these two 
quantities are fixed, scattering states can be completely 
determined by fixing in addition the scattering length of 
the 3 Si phase, and then imposing orthogonality to the 
deuteron state. This is equivalent to renormalizing with 
three counter-terms in momentum space [43| . Of course, 
more counter-terms could be considered if orthogonality 
is given up by an energy dependent boundary condition, 
as done above for the 1 S'o channel. 

In the following two subsections we will consider the 
description of the deuteron bound state and the scatter- 
ing states and discuss our results. 



B. The Deuteron 

In the case of negative energy we consider Eq. 
with 



- 7 2 = -AfB rf , 



(45) 



with 7 the deuteron wave number and Bd the deuteron 
binding energy. We solve Eq. ([55]) together with the 
asymptotic condition at infinity 



u(r) 
w{r) 



A D e-^ r 1 + 



(46) 



where As and Ajj are the s- and d-wave normalization 
factors. The asymptotic D/S ratio parameter 77 is defined 
as ?/ = A D /A S . 
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Set 


7 (fm" 1 ) 


V 


As (fm- 1/2 ) 


r m (fm) 


Qd(fm 2 ) 


Pd (%) 


(r- 1 ) 


<r- 2 } 


LO 


Input 


0.02633 


0.8681(1) 


1.9351(5) 


0.2762(1) 


7.31(1) 


0.486 (1) 


0.434(3) 


NLO-A 


Unbound 


- 


- 


- 


- 


- 


- 


- 


NLO-A (fiA = 1.34) 


Input 


Input 


0.884(3) 


1.963(7) 


0.274(9) 


5.9(4) 


0.446(10) 


0.29(2) 


NLO-A [h A = 1-05) 


Input 


Input 


0.84(4) 


1.86(8) 


0.24(3) 


12(5) 


0.62(15) 


0.8(4) 


N2LO-A (Set IV) 


Input 


Input 


0.884(4) 


1.967(6) 


0.276(3) 


8(1) 


0.447(5) 


0.284(8) 


N2LO-A (ttN) 


Input 


Input 


0.896(2) 


1.990(3) 


0.282(5) 


6.1(8) 


0.4287(13) 


0.253(2) 


N2LO-A(Fit 1) 


Input 


Input 


0.892(2) 


1.980(4) 


0.279(5) 


5.9(9) 


0.4336(15) 


0.262(3) 


N2LO-A(Fit 2) 


Input 


Input 


0.890(2) 


1.975(3) 


0.278(5) 


5.8(9) 


0.4470(15) 


0.268(2) 


Nijmll 


0.231605 


0.02521 


0.8845 


1.9675 


0.2707 


5.635 


0.4502 


0.2868 


Reid93 


0.231605 


0.02514 


0.8845 


1.9686 


0.2703 


5.699 


0.4515 


0.2924 


Exp. 


0.231605 


0.0256(4) 


0.8838(4) 


1.971(5) 


0.2860(15) 









TABLE IV: Deuteron properties for the OPE and TPE potentials. The computation is made by fixing 7 and 77 to their 
experimental values. The errors quoted in both TPE computations reflect the uncertainty in the non-potential parameters 7, 
77 and «o only. For the OPE (LO) we take q^nn = 13.1(1). For the LEC's in the TPE calculation, by set IV we refer to the 
determination of Ref [3 (see main text). For the A case we use Fits 1 and 2 of Ref. [52l ]. Fit 1 involves the SU(4) quark-model 
relation Yia = 3(?a/(2\/2) = 1.34 for <?a = 1.26. Fit 2 takes hA = 1.05. The experimental values are taken from the following 
references: 77 from [§3j], As from [6||], r m from [65| and Qd from (see also Ref. [67J for a brief review). 



In order to obtain the regularized wave functions, we 
fix 7 and 77 to their experimental values (see table [TV]) , 
and obtain u(r) and w(r) by integrating Eq. (|36p from 
r — ► 00 to r = with (|46[) as boundary conditions. As 
can be later determined from 

/>oo 

/ dr [u{rf + w(r) 2 ] = 1, (47) 
Jo 

i.e., from demanding the deuteron normalization to be 
equal to unity. The renormalized deuteron wave func- 
tions are depicted in Fig.rjat LO, NLO-A and N2LO-A, 
and compared to the Nijmegen II wave functions [6J. The 
asymptotic normalization u — > e~ ir has been adopted 
and the asymptotic D/S ratio is taken to be 77 = 0.0256. 
One can see in Fig. [2] the appearance of an increasing 
number of oscillations in the wave function when the 
radius approaches zero. They have no appreciable ef- 
fect on the physics of the deuteron, as they happen at 
very short distances. Therefore they have a very small 
effect on the computation of deuteron observables and 
cannot be resolved by external probes at the energies for 
which the effective theory description is valid. This later 
point is shown explicitly in Ref. [43| for elastic electron- 
deuteron scattering. Perhaps the most remarkable aspect 
of the present calculation is the convergence of the pro- 
posed scheme when the A-resonance is considered, which 
of course implies that the deuteron is bound at all com- 
puted orders of the potential. This is in contrast with the 
delta-less theory, where at NLO-4& the deuteron became 
unbound (as discussed extensively in Ref. (30j). 

Here, we also compute the matter radius, which reads, 

lr 2 \ 1 f°° 

rl = V = 7 / r >M 2 + w <r) 2 ) dr > ( 48 ) 



the quadrupole moment (without meson exchange cur- 



rents) 

1 1"°° 

Qd= 2tt] r2w (r)(2^2u{r) - w(r))dr , (49) 
the £>-state probability 

/>oo 

Pd= w(r) 2 dr, (50) 
Jo 

and the inverse moments of the radius 

/>oo 

{r~ n )= r- n (u(r) 2 + w(r) 2 )dr , (51) 
Jo 

which, as is well known, appear in the multiple scatter- 
ing expansion of the 7r-deuteron scattering length. Some 
results for the inverse moments in A-less effective field 
theory can be found in Refs. EH . In Tab. II VI we show 
our results in a variety of situations. In general, the re- 
sults for NLO-A and N2LO-A are in agreement with the 
experimental value of the deuteron observables, with the 
exception of the quadrupole moment, which presents a 
discrepancy of 0.01 fm 2 . The reason for it lies in meson 
exchange current contributions to the quadru pole mo- 
ment, which were estimated for ChPT in Ref. [S3], and 
are of the order of the difference of our results with re- 
spect to the experimental value. When comparing the 
A-full computations with the A-less ones, one can no- 
tice, in the first place, that a renormalized result exists 
for NLO-A, unlike the NLO-^ case (as was explained 
in more detail in Ref. 30]). This supports the better 
convergence properties of effective theory when includ- 
ing the A degree of freedom. It can also be noticed that 
the D-wave probability Pd, although not an observable, 
it is better described in N2LO-A than in N2LO-^, when 
compared with the results coming from the phenomeno- 
logical potentials [s3] . The consequence is that the mag- 
netic moment of the deuteron would be better described 
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in the A-full theory than in the A-less theory, although 
no actual computation has been made on that respect in 
the present paper. The reason for that is that in a non- 
relativistic framework, the deuteron magnetic moment 
depends solely on Prj [2Q( ■ Before comparing the results 
for other observables, one important comment must be 
made: one should only directly compare the N2LO-A 
results with the N2LO-A (irN) ones. Comparison with 
N2LO-A set IV should be done with care, as for this case 
the LECs are fitted to reproduce NN data [11]. This is 
why they look slightly better than the other results of 
Table IIVI Therefore, the comparison of N2LO-A with 
N2LO-A (irN) results implies in particular that the in- 
clusion of A enhances the compatibility between ttN and 
NN scattering, as one would expect within the EFT phi- 
losophy. 



C. Phase Shifts 

Finally, in the case of positive energy we consider 
Eq. (J36J) with 



M 



(52) 



where k is the corresponding cm. momentum. We solve 
Eq. (|36|) for the two linear independent scattering states, 
which are usually labelled as the a and (3 states. They 
are defined by the asymptotic normalization 



Uk,u(r) 

Wk,a{r) 

Wk,p(r) 



«os e ^ jo (kr) cot Si - y (kr) 
sin e (j 2 (kr) cot Si - y 2 (kr) 



-sin e 



cose 



(jo(kr) cot S 2 - ya(kr) 
(j 2 (kr) cot S 2 - y 2 (kr) 



(53) 



(54) 



where ji (x) — xji (x) and yi (x) — xyi (x) are the reduced 
spherical Bessel functions and 8\ and 8 2 are the eigen- 
phases in the 3 Si and 3 Di channels; e is the mixing angle 
Ei . The orthogonality constraints between the deuteron 
and scattering states generate the following boundary 
conditions 



u l u 'k,a ~ Uy11>k,a + W-fW k ,a ~ ^Wfc, Q 
u l u 'k,l3 - u y u k,f3 + W^w' kJ3 - w' 7 Uk,f3 



0. 



= 0. 



(55) 



The use of the superposition principle for the a and (3 
scattering states, plus the deuteron wave functions, al- 
lows to deduce the corresponding 3 Si — 3 Di eigen phase- 
shifts. The results are depicted in Fig. [3] at LO, NLO-A 



and N2LO-A compared to the Nijmegen II potential re- 
sults |6[. We observe a clear improvement in the thresh- 
old region, and quite remarkably for the Ei mixing phase. 
However, there is no improvement in the 3 Di phase-shift. 



COMPARISON WITH THE DELTA-LESS 
THEORY 



In our previous work [301 ] , the renormalization of Delta- 
less theory was analyzed. As discussed above, it was 
found that at NLO-A there was no deuteron bound state 
if the cut-off was removed. The reason was due to the 
short distance ~ g\/(f^r 5 ) repulsive singular character 
of the potential in the 3 Si — 3 Di channel. However, at 
N2LO-A three counter-terms where needed due to the 
short distance attractive singular character of the po- 
tential. In fact the agreement with more sophisticated 
calculations [24[ was remarkable. More recently, the 
quality of these chiral wave functions has been tested in 
electron-deuteron scattering [431 ] using LO currents, with 
an amazingly good agreement up to momentum transfer 
of q = lGeV. It is interesting to compare the results of 
the Delta-less theory [3(| with the ones found here af- 
ter inclusion of the A in the potential. The first aspect 
to note is that the NLO-A potential has an attractive 
~ g\l (A/^r 6 ) short distance singularity, a feature kept 
at the N2LO-A, however with different scales. 

In Fig.[4]the np 1 5o singlet phase shifts are depicted as 
a function of the cm. momentum for the N2LO-A and 
N2LO-A compared to the Nijmegen II potential phase 
shifts 6|. As we see, with only one counter-term, i.e. fix- 
ing the scattering length ao, the result is slightly wors- 
ened when the A is included. This is consistent with the 
change in the effective range reported in Table IIIII Of 
course, if ro is fixed as an additional renormalization con- 
dition, there is an improvement in the low energy region 
but the difference between including or not A degrees of 
freedom is hardly visible. 

The situation for the deuteron wave functions is 
slightly different. As we see in Fig.O the present N2LO- 
A deuteron wave functions resemble slightly better the 
Nijmegen II ones 0, In particular, the D-wave be- 
comes smaller when one goes from our previous TPE 
ones (30j . 

Finally, in Fig. [6] we show the np eigen phase shifts 
in the 3 S'i — 3 Di channel as a function of the cm. mo- 
mentum for the N2LO-A and_N2LO-A compared to the 
Nijmegen II potential results [6| ■ The noticeable improve- 
ment in the Ei phase is in agreement with the smaller 
D-wave deuteron wave function. 
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to the Nijmegen II potential phase shifts [f|. 
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VI. SCHEME DEPENDENCE AND CUT-OFFS 
A. The relevant scales 

As we see our scheme including the A does not repro- 
duce the 1 Sq phase shift for momenta larger than the pion 
mass even if the effective range is adjusted to its experi- 
mental value. This trend has also been observed in other 
renormalized calculations including TPE effects and 1 /M 
corrections in a Heavy Baryon formalism [3fj| . using a 
relativistic potential [71| or including N3LO contribu- 
tions [42| . Within the chiral approach to NN interactions 
the main candidates for explaining such a disagreement 
would be 37r-exchange or relativistic corrections. 

It is interesting at this point to ask which are the rel- 
evant scales which build the full strength of the phase 



shifts. As already found in [29|, |30|, |3l|, |42j, lZl| this is 
about 0.4 — 0.5 fm for Delta-less calculations. For illus- 
tration purposes, the short distance cut-off dependence is 
shown in Fig.[7]for the effective range in the case of renor- 
malization with one counter-term in the Heavy Baryon 
Delta-less theory to LO, NLO-4& and N2LO-^ where also 
the Nijmegen II potential [(| is considered. As we see, and 
regardless on the full renormalized effective range value, 
the Nijmegen II saturating scale is in between NLO-4& 
and N2LO-A., a not unreasonable result. Likewise the 
N2LO-A effective range displays a similar approach to 
the renormalized result. 

These scales are comparable to the nucleon size but 
also to the range where 37r exchange starts contributing 
since it behaves as e - 3m ^ r at long distances. According 
to the results of Kaiser (72l . 1731 . Lz3 (and in the absence 
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FIG. 6: The np spin triplet (eigen) phase shifts for the total angular momentum j = 1 as a function of the cm. momentum 
for the N2LO-A and N2LO-^& potentials compared to the Nijmegen II potential results [!, H|]. 



of A) the potential is attractive, and using dimensional 
estimates for the 37r-exchange potential, it behaves as 
g\/{f^r 7 ) at short distances indicating a stronger singu- 
larity and therefore a stronger short distance u(r) ~ r* 
suppression as well. Thus, we expect that including these 
effects would only slightly worsen the results by reducing 
the phase-shift. 

Another interesting observation is that by taking larger 
coordinate space cut-offs the phase shift is not improved, 
as in our regularization procedure it converges from be- 
low, so in this case the best possible cut-off is r c — > 9 . 



9 Improving or not the phase shifts by using a finite cut-off while 
keeping the same renormalization conditions is, in fact, a regular- 
ization scheme dependent feature, and hence a further good rea- 
son to remove the cut-off. This choice however does not resolve 
the representation dependence of the potential on the choice of 
the pion field. Therefore, the cut-off less solution will obviously 



In this context, one can identify two kind of finite cut-off 
sources of errors: the first one is related with the explicit 
form of the regulator employed in the computations, and 
the second one with the actual size of the finite cut-off. 
This second source of errors is the only one which is usu- 
ally assessed in most effective field theory works, while 
the role played by the regulator is commonly ignored, 
probably resulting in an underestimation of the errors. 
While in our present renormalization scheme the most 
sensible thing to do in the 1 S'o singlet channel is to com- 
pletely remove the cut-off, this may be not be the case 
in other regularization schemes. Although this could be 
in fact considered as a good motivation for keeping a fi- 
nite cut-off or introducing form factors, we think that 
this kind of procedure is difficult to justify from the EFT 



depend on the previous choice. 
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viewpoint. 

The influence of relativistic effects on the results is less 
obvious, since they generate either energy dependence or 
non-localities, but experience in renormalizing relativis- 
tic potentials indicates that they are not crucial, at least 
in the 1 Sq channel [7l|, and are saturated by scales larger 
than the nucleon Compton wavelength, i.e. well above a 
possible influence from the NN channel. 

Thus, we see that in all these calculations including 
TPE effects the saturating scales are of the order of 
0.5 fm, regardless on detailed issues, in particular includ- 
ing or not the A. This is comparable to other scales, 
2/m p ~ 0.51 fm, which may equally represent vector 
meson exchange or nucleon size effects. Therefore the 
discrepancy between renormalized phase shifts and phe- 
nomenological ones should be considered a real one, in 
the sense that finite size effects may be important physi- 
cally. Moreover, it should be kept in mind that once the 
short distance cut-off becomes smaller than the nucleon 
size, the effect of the counter-terms is not enough to re- 
produce intermediate energy phase shifts, meaning that 
they cannot mimic the finite nucleon size or other short- 
distance physical effects which may be responsible of the 
observed discrepancy. 



B. Remarks on Spectral Regularization 

The calculation of the NN potentials can advanta- 
geously be carried out by using the method of disper- 
sion relations [HI, [lj, This has motivated the use 
of the so-called spectral regularization [26|, [27J, [28|, [52| in 
NN calculations where finite momentum space cut-offs 
have been implemented. Remarkably, a good fit to the 
1 S'o phase was achieved. This is in contrast to the N3LO 
computation carried out in Ref. [42] with one countert- 
erm, where a systematic underestimation of the data was 
found. In our view, the relevant issue is to disentangle 
the cut-off artifacts from clearly attributable physical ef- 
fects. In this section we analyze the issue of cutting-off 
the potential. 

The two pion exchange potential satisfies a dispersion 
relation based on the NN — > 2ir amplitude. One has the 
representations 



and r <C A 1 , the potentials behave as 



Vc(r,A) = / dnnp c <ji)e ^ , 

-" ' Jim 

r A 



A 



V s (r,A) 



V T (r,A) 



1 



6n 2 r 



g 7r 2 r 3 



2m 



dpi* {p 2 Pt{p) -Sps(fJ-)) e Air , 



dpp 3 {3 + 3pr + p 2 r 2 )p T (p)e'^ r 



2m 



(56) 



with Pi(p) — ImVi{ip) and similar relations for the Wj 
potentials. It is straightforward to see that for small A 



V c (r,A) 
V s (r,A) 
V T (r,A) 



2ir 2 r 



dpppc{p) , 



g 7r 2 r 3 



2m 
A 



dpp (p 2 p T {p) - 3ps(jJt)) , 



dpp 3 3p T (p) . 



2m 



(57) 



Thus, although Vc, Wc, Vs and Ws become regular, the 
tensor contributions Vt and Wt remain singular, despite 
higher energy states being cut off. That means that reg- 
ularization of the scattering problem in triplet channels 
is mandatory to obtain well defined results. 

Note that this short distance behaviour is quite dif- 
ferent from the one obtained for r > when A — > 00 
(see appendix EJ. In Fig. [5] the ratio, V(r, A)/V(r), of 
1 Sq potentials as they enter in the phase shift calculation 
with and without spectral regularization is depicted. We 
note that there is a sizeable distortion at 2-3 fm's with 
A = 700 MeV, decreasing the strength of the interac- 
tion and hence providing effectively a repulsion in the 
singlet potential. The effects of the spectral regulariza- 
tion on the potential can be seen in Fig. [5] for NLO-A 
and N2LO-A both with one counter-term as well as with 
two counter-terms. As we see, it is possible to describe 
the data successfully for suitable values of the spectral 
cut-off A ~ 700 McV, in agreement with the findings of 
Ref. 52j. It is noteworthy that this happens regardless 



on the usage of one or two counter-terms, reinforcing the 
conclusion that the agreement is mainly due to the dis- 
tortion in the potential. In regard to the behaviour with 
respect to spectral regularization, the trends presented 
for the A-full theory are also reproduced in the A-less 
scheme. We remind that an accurate description of the 
data was achieved in Ref. [28[ at N3LO-A with a spec- 
tral regularized potential with finite momentum space 
cut-offs and 4 counter-terms in the 1 S'o channel (the full 
N3LO computation has 24 counter-terms) . On the other 
hand the data were not described for large values of p in 
the N3LO-A and (s pec trally) unregularized renormalized 
calculation of Ref. [43] • 

As before, it is interesting to analyze the relevant scales 
building up the effective range in the case of renormal- 
ization with one counter-term. In Fig. [3 we show the 
effective range ro as a function of the short distance cut- 
off r c for several values of the spectral cut-off A. Clearly, 
the stability region in r c is shifted towards lower values 
as the spectral cut-off is decreased. This is consistent 
with the large distortion of the potential at intermediate 
scales. 

We have also analyzed the impact of the spectral reg- 
ularization on the already successful description of the 
deuteron presented in previous papers and the present 
one. The results do not change noticeably after (spec- 
trally) regularizing the potential: they lie between the 
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FIG. 7: 1 5'o effective range (in fm) ro(r c ) — 2 f / °°(1 — r/ao) 2 dr — J r °° uq drj , with ao = — 23.74fm as a function of the short 
distance cut-off r c (in fm). We compare the (Heavy Baryon) Delta-less theory at LO, NLO and N2LO and the Nijmegen II 
potential [f| (left panel) with the spectral-regularized N2LO-A potential for different values of the spectral cut-off A (right 
panel) . 



5 0.6 
> 

is 0.4 



1 .J*. " 




' ' / 
• / / 


' / / 




/ /// V sr /V full 




: i / / 




' i / 




// / 




, i 




' 1 / 




;// 600 MeV 

V / 700 MeV ■ 






,'// 800 MeV 


\ jy 2000 MeV 



1 2 3 4 5 



r[fm] 

FIG. 8: V(r, A)/V(r)-ratio of 1 So potentials as a function of 
the distance r (in fm) for the spectral-regularized N2LO-A 
potential with different values of the spectral cut-off A as a 
function of the short distance cut-off r c (in fm). 



results obtained with the OPE potential and with NLO- 
and N2LO-A. This is compatible with the weakening 
of these contributions to the chiral potential. Although 
the final results are rather simple to summarize, some re- 
marks should be added on the renormalization procedure 
for the deuteron channel when spectral regularization is 
applied. The modified NLO- and N2LO-A chiral poten- 
tials have now an attractive and repulsive eigen-channel, 
instead two attractive ones, which lead us to the following 
alternative: either we remove the cutoff, in which case we 
can only fix the binding energy and obtain the D-to-S ra- 
tio rj as a prediction, or we keep a finite cutoff, and then 



fix both the binding energy and rj. Provided that the 
finite cut-off is sensible, about r c ~ 0.5 fm (i.e. the sat- 
uration scale) for a spectral cutoff A ~ 600 — 800 MeV, 
both procedures give equivalent results. The previous 
discussion on the effects of spectral cut-off agree with 
the remarks presented in Ref. [43| regarding the effects 
of this regularization of the potential for the deuteron 
form factors. For completeness we show the results in 
Table E0 

In summary, reducing the strength of the potential at 
short distances by means of the spectral regularization is 
phenomenologically preferred in the singlet channel and 
innocuous in the triplet channel, but distorts largely the 
chiral potential in a region much larger than the nucleon 
size. Moreover, we find that within our scheme the agree- 
ment with the data is achieved regardless of the addi- 
tional counter-terms invoked by Weinberg's counting, be- 
ing in fact superfluous after renormalization. This said, 
although the spectral regularization improves over the 
standard finite cut-off approaches and seems phenomeno- 
logically favoured, the inclusion of finite size effects in a 
model independent manner would certainly be very use- 
ful. 



VII. CONCLUSIONS 

Since more than fifteen years there has been a growing 
interest in pursuing an EFT description of NN scattering 
below pion production threshold where the spontaneous 
breakdown of chiral symmetry in QCD is manifestly ex- 
ploited. One of the reasons which have greatly hindered 
the EFT developments within the NN context has been 
the lack of a credible power counting for the potential 
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Set 


7 (fm" 1 ) 


V 


As (fm~ 1/2 ) 


r m (fm) 


Qd(fm 2 ) 


Pd (%) 


(r- 1 ) 


(r- 2 ) 


LO 


Input 


0.02633 


0.8681(1) 


1.9351(5) 


0.2762(1) 


7.31(1) 


0.486 (1) 


0.434(3) 


NLO-A 


Unbound 


- 


- 


- 


- 


- 


- 


- 


A = 700 MeV 


Input 


0.02669 


0.851(7) 


1.894(15) 


0.270(4) 


8.9(1.0) 


0.6(2) 


1.8(1.5) 


r c = 0.5 fm 


Input 


Input 


0.83(2) 


1.86(6) 


0.24(2) 


8(3) 


0.56(7) 


0.51(15) 


NLO-A (h A = 1.34) 


Input 


Input 


0.884(3) 


1.963(7) 


0.274(9) 


5.9(4) 


0.446(10) 


0.29(2) 


A = 700 MeV 


Input 


0.02637 


0.8720(10) 


1.938(2) 


0.2781(6) 


7.31(12) 


0.48(2) 


0.5(2) 


r c = 0.5 fm 


Input 


Input 


0.867(13) 


1.93(3) 


0.263(14) 


6.6(1.2) 


0.48(4) 


0.35(8) 


NLO-A (h A = 1.05) 


Input 


Input 


0.84(4) 


1.86(8) 


0.24(3) 


12(5) 


0.62(15) 


0.8(4) 


A = 700 MeV 


Input 


0.02651 


0.864(2) 


1.922(5) 


0.2755(14) 


7.8(3) 


0.52(5) 


0.7(4) 


r c = 0.5 fm 


Input 


Input 


0.854(17) 


1.90(4) 


0.26(2) 


7(2) 


0.51(5) 


0.41(11) 


N2LO-A (Set IV) 


Input 


Input 


0.884(4) 


1.967(6) 


0.276(3) 


8(1) 


0.447(5) 


0.284(8) 


A = 700 MeV 


Input 


0.02504 


0.877(2) 


1.946(5) 


0.264(2) 


5.6(7) 


0.48(2) 


0.42(6) 


r c = 0.5 fm 


Input 


Input 


0.875(8) 


1.942(15) 


0.270(2) 


8(2) 


0.46(2) 


0.27(6) 


N2LO-A (ttN) 


Input 


Input 


0.896(2) 


1.990(3) 


0.282(5) 


6.1(8) 


0.4287(13) 


0.253(2) 


A = 700 MeV 


Input 


0.02590 


0.8801(2) 


1.9540(5) 


0.27659(14) 


6.580(6) 


0.463(4) 


0.35(3) 


r c = 0.5 fm 


Input 


Input 


0.883(4) 


1.956(6) 


0.274(9) 


5.9(8) 


0.447(13) 


0.28(2) 


N2LO-A(Fit 1) 


Input 


Input 


0.892(2) 


1.980(4) 


0.279(5) 


5.9(9) 


0.4336(15) 


0.262(3) 


A = 700 MeV 


Input 


0.02603 


0.8783(3) 


1.9507(6) 


0.2774(2) 


6.77(3) 


0.465(4) 


0.36(4) 


r c = 0.5 fm 


Input 


Input 


0.880(6) 


1.954(14) 


0.272(10) 


5.9(6) 


0.45(2) 


0.29(4) 


N2LO-A(Fit 2) 


Input 


Input 


0.890(2) 


1.975(3) 


0.278(5) 


5.8(9) 


0.4470(15) 


0.268(2) 


A = 700 MeV 


Input 


0.02606 


0.8769(4) 


1.9479(7) 


0.2770(2) 


6.83(3) 


0.467(5) 


0.37(4) 


r c = 0.5 fm 


Input 


Input 


0.878(6) 


1.950(15) 


0.271(10) 


5.9(6) 


0.46(2) 


0.30(4) 


Nijmll 


0.231605 


0.02521 


0.8845 


1.9675 


5.635 


0.2707 


0.4502 


0.2868 


Reid93 


0.231605 


0.02514 


0.8845 


1.9686 


0.2703 


5.699 


0.4515 


0.2924 


Exp. 


0.231605 


0.0256(4) 


0.8838(4) 


1.971(5) 


0.2860(15) 









TABLE V: Deuteron properties for the OPE and TPE potentials with spectral regularization. The computation is made by 
fixing 7 and rj to their experimental values or by fixing 7 and predicting rj depending on the singularity structure of the 
potential and the coordinate space cut-off. For each case we present three computations: (i) the complete computation, i.e. 
the results obtained when the cut-off is completely removed and there is no spectral cut-off, (ii) the spectrally regularized 
computation with r c — 0.1 fm, and (iii) the spectrally regularized computation with r c = 0.5 fm. The errors quoted in the 
(spectrally) unregularized TPE computations and in the spectrally regularized potential with a finite cut-off r c — 0.5 fm reflect 
the uncertainty in the non-potential parameters 7 and rj. In the spectrally regularized potentials, the errors represent the cut-off 
dependence of the results for r c ranging from 0.1 to 0.2 fm. For the OPE (LO) we take g^NN = 13.1(1). We take set IV [24[ 
for the LEC's in the TPE calculation. For the A case we use Fits 1 and 2 of Ref. [52| . Fit 1 involves the SU(4) quark-model 
relation h A = 3g A /(2V2) = 1.34 for g A = 1.26. Fit 2 takes h A = 1.05. 



which at the same time complies to short distance in- 
sensitivity. Given the tight constraints under which this 
might actually happen, it has not been obvious which 
particular form of the chiral expansion indeed embod- 
ies these desirable properties. In the present work we 
have shown how the inclusion of A degrees of freedom in 
the potential not only complies to a phenomenologically 
well founded motivation, but also provides the requested 
short distance insensitivity of the central phases and the 
deuteron properties after the necessary renormalization is 
carried out. This improves the previous situation without 
explicit A degrees of freedom: while at and N2LO- 

4& the existence of a deuteron was compatible with renor- 
malizability, at that was not the case. This has 

raised reasonable doubts on the suitability and useful- 



ness of non-perturbative renormalization per se to chiral 
potentials. A very rewarding aspect of the present in- 
vestigation is the existence of a deuteron bound state at 
LO, NLO-A and N2LO-A. Of course, a proof of consis- 
tency to all orders, including relativistic, spin-orbit, three 
pion exchange corrections, etc., remains to be done. This 
is so because, although the power law behaviour of the 
chiral potential at short distances can be trivially fixed 
by dimensional arguments and power counting, the de- 
termination of the attractive-repulsive character of the 
potential can so far only be fixed by actual calculations. 
The found deuteron properties seem to obey a converg- 
ing pattern and in the E± phase a clear improvement is 
observed at N2LO-A. Moreover, our N2LO-A results re- 
semble much those of N2LO-^ after renormalization for 
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FIG. 9: np 1 So renormalized phase shifts for the spectral-regularized N2LO-A potential for different values of the spectral cut- 
off A as a function of the cm. momentum (in MeV) compared to the Nijmegen II potential results 0] with one counter-term 
(left panel) and with two counter-terms (right panel), corresponding this last choice to the standard Weinberg counting. 



reasonable parameter values describing the ttN reaction 
close to threshold. This suggests that despite the previ- 
ous inconsistency in NLO-4&, the N2LO-^ deuteron wave 
functions can be used for practical purposes, despite the 
theoretically unpleasant and disturbing "jumping" of the 
calculation. In addition, it should be reminded 
that within these approximations the ttN threshold prop- 
erties can be fitted, with the sole exception of b$ + , and 
thus the overall consistency between the ttN and NN 
sectors is almost satisfied. Therefore, the inclusion of the 
A resonance complies with the original EFT motivation 
of describing simultaneously nN and NN scattering at 
low energies. 

We have also found that at the level of approximations 
involved in the present and previous renormalized calcu- 
lations, the 1 Sq phase shift is not entirely reproduced for 
cm. momenta larger than the pion mass if we insist on a 
reasonable ttN physics. Several effects might be respon- 
sible for this persistent discrepancy. We have discussed 
those on the light of the relevant scales, r > r c = 0.5 fm, 
which practically provide the total contribution to ob- 
servables. We have further discussed parameterizations 
of the NN force based on a spectral regularization of the 
NN potential with a cut-off of A = 700 MeV. We have 
shown that agreement to data is achieved mainly due to a 
large distortion of the potential at 2-3 fm and regardless 
on additional inclusion of counter-terms. These length 
scales are much larger than the expected finite nucleon 
size or vector meson exchange effects, casting doubts on 
the usefulness of spectral regularization and suggesting 
the need for a more controllable and better founded de- 
scription of the missing short distance physics in the 
channel. 

In the present paper we have restricted to central 
phases and the deuteron, but the calculation of higher 
partial waves is also of interest, as well as the inclusion 
of Coulomb effects in pp scattering. In all cases, the neg- 



ative definite character of the potential at short distances 
guarantees the existence of convergent results. Finally, 
the present results can have some impact on calculations 
including deuteron properties such as deuteron form fac- 
tors, pion-deuteron scattering, and further low energy 
matrix elements of electroweak deuteron reactions where 
short distance insensitivity and chiral symmetry are both 
expected to play a significant role. 
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APPENDIX A: SHORT DISTANCE EXPANSION 
OF THE POTENTIALS 

Using the spectral representation the short distance 
expansion of the potentials can be done by expanding 
spectral functions for fi»m, with \ir fixed. This way 



19 



one obtains 



and at N2LO-A by 
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